function [log_like,log_like_deriv]=find_log_like(theta,mat1)

	eps=1e-6;
	KK=length(theta);
	mat2=mat1*theta;
	log_like=sum(log(mat2 + eps),1);
	log_like_deriv=sum((1./(mat2 + eps)*ones(1,KK)).*(mat1),1);

	%

